Decoherence in weak localization II: Bethe-Salpeter calculation of Cooperon 
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This is the second in a series of two papers (I and II) on the problem of decoherence in weak 
localization. In paper I, we discussed how the Pauli principle could be incorporated into an influence 
functional approach for calculating the Cooperon propagator and the magnetoconductivity. In the 
present paper II, we check and confirm the results so obtained by diagrammatically setting up a 
Bethe-Salpeter equation for the Cooperon, which includes self-energy and vertex terms on an equal 
footing and is free from both infrared and ultraviolet divergencies. We then approximately solve 
this Bethe-Salpeter equation by the Ansatz C{t) = & {t)e~^''*\ where the decay function F{t) 
determines the decoherence rate. We show that in order to obtain a divergence-free expression 
for the decay function F{t), it is sufficient to calculate C^{t), the Cooperon in the position-time 
representation to first order in the interaction. Paper II is independent of paper I and can be read 
without detailed knowledge of the latter. 



I. INTRODUCTION 

This is the second in a series of two papers (I and II), 
in which we revisit the problem of decoherence in weak 
localization, using both an influence functional approach 
(paper I) and a Bethe-Salpeter equation for the Cooperon 
(paper II) to calculate the magnetoconductivity. The ba- 
sic challenge is to calculate the interference between two 
time-reversed trajectories of an electron travelling diffu- 
sively in a Fermi sea and coupled to a noisy quantum 
environment, while taking proper account of the Pauli 
principle. In paper li, we discussed how this could be 
done using an influence functional approach by dress- 
ing the spectrum of the noise field by "Pauli factors" 
[see Eq. (II66|): throughout, "I" will indicate formulas 
from paper I]. Moreover, within the influence functional 
scheme we concluded that a divergence-free calculation 
of the decoherence rate can be obtained by expressing 
the Cooperon in the position-time representation as 



C(0,t) ~(7°(0,^)e-^(*^ F{t) 



(70(0, i) 



(1) 



where C^(0,t) is the first-order term in an expansion of 
the full Cooperon C(0,t) in powers of the interaction. 
[In the present paper, this statement will be made more 
precise: when reexponentiating, a part of C^(0, i) has 
to be omitted that can be determined, in a self-energy- 
only calculation, to contribute only to the prefactor of 
the Cooperon, see Sec. Ill 

These conclusions of paper I rested entirely on the in- 
fluence functional approach, and, in the discussion of the 
Pauli principle, relied on heuristic arguments. Though 
these are in accord with results derived elsewhere^iMi^ 
(as shown in Paper 1, Section lVlI|) . it is desirable to com- 
pare the approximations used and the results obtained so 
far against a treatment relying purely on diagrammatic 
perturbation theory, the framework within which most of 



our understanding of disordered metals to date has been 
obtained. 

In the present paper II, we check and confirm the re- 
sults mentioned above by diagrammatically setting up 
a Bethe-Salpeter equation for the Cooperon using stan- 
dard Keldysh diagrammatic perturbation theory (using 
conventions summarized in^), which includes self-energy 
and vertex terms on an equal footing and is free from 
both infrared and ultraviolet divergencies. We then show 
that this equation can be solved (approximately, but with 
exponential accuracy) with an Ansatz that is precisely of 
the form Eq. Q, and that the function F{t) so obtained 
agrees with the form derived in paper I [Eq. (I l()5|l ]. 

The usual diagrammatic calculation of the Cooperon 
starts from a Dyson equation for a "self-energy-diagrams- 
only" version of the Cooperon, 



c::;(c.) = c^(c.) [i-hi]^4V)c^fH] 



(2) 



Here the Cooperon self-energy S|®^^ includes only self- 
energy diagrams, in which interaction lines connect only 
forward to forward or backward to backward electron 
propagators; for these diagrams, the frequency labels 
along both the forward and the backward propagators 
are conserved separately, which is why the Dyson equa- 
tion is a simple algebraic equation for C|°J^(aj). However, 
the Cooperon self-energy S|°J^(w) turns out to be infrared 
divergent in the quasi 2- and 1-dimensional cases. This 
problem is usually cured by inserting an infrared cutoff 
by hand (as reviewed in Section III CI below) . The results 
so obtained are qualitatively correct but, due to the ad 
hoc treatment of the cutoff, not very accurate quantita- 
tively [e.g. in the first hnc of Eq. ijEJ for "'(<) below 
the exponent is correct, but the prefactor is wrong by 
roughly a factor of 2 compared to Eq. (l E|l ]. 

Our goal in the present paper is to obtain, starting 
from a diagrammatic equation, results free from any cut- 
offs, infrared or ultraviolet, that have to be inserted by 



2 



hand - the theory "should take care of its divergencies 
itself" . This can be achieved if the Cooperon self-energy 
is taken to include vertex diagrams, in which interac- 
tion lines connect forward and backward electron prop- 
agators. Since this brings about frequency transfers be- 
tween the forward and the backward propagators, their 
frequency labels are no longer conserved separately. As a 
consequence, it becomes necessary to study a more com- 
plicated version of the Cooperon, (ili, 172), labeled by 
three frequencies, and governed not by a simple algebraic 
Dyson equation, but by a nonlinear integral equation, 
which we shall refer to as "Bethe-Salpeter equation" . 

Finding an exact solution to the Bethe-Salpeter equa- 
tion seems to be an intractable problem, which we shall 
not attempt to attack. Instead, we shall transcribe the 
Bethe-Salpeter equation from the momentum-frequency 
to the position-time domain, in which it is easier to make 
an informed guess for the expected behavior of the solu- 
tion. Using the intuition developed in paper I within the 
influence-functional approach [summarized in Eq. Q of 
the present paper] , we shall make an exponential Ansatz 
for C^(ri2, ii, i2), the Cooperon in the position-time do- 
main. We shall show that this Ansatz solves the Bethe- 
Salpeter equation with exponential accuracy, in the sense 
that improving the Ansatz would require terms to be 
added to the exponent that are parametrically smaller 
(in powers oil/ g, g being the dimensionless conductance) 
than the leading term in the exponent. 



II. SETTING UP BETHE-SALPETER 
EQUATION FOR COOPERON 
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Figure 1: (a) Diagram for the weak localization correction to 
the AC conductivity, 5a^^(uJo) [Eqs. Q]. In contrast to the 
so-called "interaction corrections" to the conductivity, each 
current vertex is attached to both a retarded and advanced 
electron line, (b) Diagrammatic definition of full Cooperon 
Cq(VLi,VL2), and schematic depiction of the Bethe-Salpeter 
equation JSJ satisfied by it; £ = + '® conserved 

average of the energies of the upper and lower lines, while fJi 
and 0,2 axe the outgoing and incoming Cooperon frequencies 
(with Q,j — Slj" — 0,~). For the structure of the Cooperon 
self-energy E*"" and details of our diagrammatic conventions, 
see Appendix 1X1 Fig. |21 



A. Various expressions for conductivity 

The diagrammatic definition of the weak localiza- 
tion contribution to the AC conductivity of a quasi d- 
dimensional disordered conductor is given by Fig. Q^a), 
which corresponds to the following expression (cf. 
App. A of Ref. or App. C of Ref. @]): 



cond 



(^o) 



(3a) 

■w,Wo-l-w), (3b) 



where 



/(£_)-/(£+) 



(4) 



denotes an average over e, with e± = e ± ^wo, and 
in the DC limit wq — *■ the weighting function re- 
duces to —/'(e), the derivative of the Fermi function 
/(e) = \l\e^l^ + 1]. (In this paper, temperature is 
measured in units of frequency, i.e. T stands for k^T/h 
throughout; likewise, although e will often be referred to 
as "excitation energy", it stands for a frequency.) 



The full Cooperon with general arguments, (fii, VI2) 
is defined diagrammatically in Fig. ^b): E is the aver- 
age of the frequencies of the upper and lower electron 
lines, while f^i and ^2 ai'e the outgoing and incoming 
Cooperon frequencies, respectively. In the absence of 
external time-dependent fields, the average energy £ is 
conserved between incoming and outgoing lines. The 
Cooperon needed for the AC conductivity in Fig.^a) has 
incoming upper and lower electron lines with energies e+ 
and E- — uj, and outgoing upper and lower electron lines 
with energies e+ — tZ; and e_, implying S7i = (wq ~ <^)j 

0. 2 = {1^0 + '^) and £ = £ ~ ^uj, as used in Eq. (|3bp . 

To make contact with the expression for the conduc- 
tivity in the position-time representation used in paper 

1, we rewrite Eq. I|3b|l as 



''cond 



dt e'"°*C'^+(0,i) 



(5a) 



C'+{0,t) = C"+(ri2 =0:ti = \t,t2 = -\t) , (5b) 



where 6*^+ (f 12; ^i, ^2) is a representation of the full 
Cooperon in an energy/position/two-time representa- 
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tion, 



(5c) 



(dq) {d 2lI))idLu) e»[9'-i2— '*i2+^*i2] 



X C 



{UJ — LU , LU + LU) , (5d) 



where ri2 = ri - r2, ti2 = - ^2, ti2 = + t2, ^ = 
i(r22 + ill) and cj = ^{fl2 — ^i)- The J dt time integral 
in Eq. H5a|l sets ui — ujq and hence e+ — ^{uj + uj) = s—}^u) 



in Eq. (|5c|) . as needed for Eq. (j3b|l . Inserting Eq. ISa] 
into (I3al) we find 



die'""* (C"'+(0,t) 



(6) 



Thus, the DC hmit 5cr^^(0) is seen to be an energy- 
averaged version of Eq. (IQ). Since our goal is to make 
contact with the results of paper I, we shall take the 
DC limit t^o ^ and e+ — s- e throughout below, (it is 
straightforward to reinstate the wo-dependence explicitly 
by replacing the parameter e by e+ in all Cooperons be- 
low). 

We choose to normalize the full Cooperons such that 
in the absence of interactions, they reduce as follows to 
their noninteracting versions {^12 = fli — ^2)'- 



Cf(fli,r!2) 
C'^(ri2;ti,t2) 



27r5(r!i2)C°(r!i) , 



no lilt 



C^(ri2,ii2) 



(7a) 
(7b) 



Here C°(ll) = {Eq - iO)-\ with Eq = Dq^ + 7//, where 
7h is a magnetic-field induced decay rate. For later ref- 
erence, we also define E^ — Dq^. 

Our strategy for determining the decoherence rate will 
be to find an approximation for C^(0, t) of the form 
To this end, we shall set up a Bethe-Salpeter equation 
for Cq{ili, ft2), transcribe it to the position-time domain 
to find a Bethe-Salpter equation for {ri2;ti,t2), and 
then solve the latter using the Ansatz 



B. Bethe-Salpeter Equation for Cooperon 



In the presence of interactions, the full Cooperon 
Cq{il.i , ^,2) satisfies a Bethe-Salpeter equation of the gen- 



eral form 

C^(fli,r!2) = C"(fli) 



2n6iQi2) 



(8) 



(dr!3)Sf_fuu(f^l,f^3)C^(fl3,fl2) 



depicted schematically in Fig. QJb) . The average energy 
£ is conserved, because no external fields are present. For 



the Cooperon self-energy E^f^[[(r2i, fia) occuring herein, 
we shall adopt the diagrammatic definition first written 
down in Ref. 7. The corresponding diagrams and equa- 
tions for Sfuii are rather unwieldy and hence have been 
relegated to Appendix IXI [see Fig.|2and Eqs. I|A2|I in Ap- 
pendix lA Ij . This very technical appendix can be skipped 
by casual readers; its contents are summarized in the next 
two paragraphs, and to be able follow the developments 
of the main text below, it should suffice to just occasion- 
ally consult the final formulas for the self-energies given 
in Eqs. (|X5|) . 

The Cooperon self-energy Sfuii is itself proportional to 
the Cooperon C, thus the Bethe-Salpeter equation (jHJ 
is non- linear in C. Solving it in full glory thus seems 
hardly feasible. Therefore, we shall henceforth consider 
only a "linearized" version thereof, obtained [in Ap- 
pendix lA 2j by replacing the full Cooperon self-energy 
Sfuu in Eq. © by a bare one, Sbarc- The latter, given 



explicitly in Eqs. (jA5|) . is obtained by making the re- 
placement (111, lis) ^ 2n6{ni ~ r23)C°(fli) for every 
occurence of the full Cooperon in Sfuii. 

A perturbative expansion of the full Cooperon C in 
powers of the interaction can readily be generated by it- 
erating Eq. (jHJ . This is done explicitly to second order in 
Appendix IX3I [see Ea. (|X6|) ]. The expansion illustrates 
two important points: firstly, due to the frequency trans- 
fers between the forward and backward propagators gen- 
erated by the vertex diagrams, the frequency arguments 
of Share increasingly become "entangled" from order to 
order in perturbation theory, i.e. they occur in increas- 
ingly complicated combinations. This makes it exceed- 
ingly difficult to directly construct an explicit solution. 
Secondly, no ultraviolet divergencies arise in perturbation 
theory, confirming the heuristic Golden Rule arguments 
of Paper I, Section (and contradicting suggestions to 
the contrary implicit in Refs. .8..9..ICk see Appendix lA 31 
for a discussion of this point). 



Recover Dyson Equation by Neglecting Vertex 
Terms 



Before attempting to solve the (linearized) Bethe- 
Salpeter equation, it is instructive to start for the mo- 
ment with a rather strong approximation, namely to sim- 
ply neglect all vertex terms (they will be reinstated later), 
thereby avoiding the abovementioned "entanglement" of 
frequencies. This reduces the Bethe-Salpeter equation 
to the more familiar Dyson equation (0) for the "self- 
energy-only" Cooperon C"'^, and will allow us to review 
some standard arguments and to recover some familiar 
and simple results. 

In the absence of vertex diagrams, the Cooperon 
self-energy T,^^^^^^{Qi,fl2) of Eq. (|A5a|) is propor- 
tional to (5(f2i2), implying the same for the Cooperon 

^^^'''°'^(fli, ri2), so that the Cooperons needed on the 
right-hand sides of Eqs. lfSc|) and (|5b|) can, respectively, 
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be written as 

7£— 77^2 ,sclf 



c 



C"^'^'="(0,t) 



(9a) 



{dq){duj) e 



(9b) 



The "single-frequency" Cooperon C| ^ (w) introduced in 



Eq. (|9a|l is the generahzation of the free, single-frequency 
Cooperon C°((jj) to the case of a Cooperon for pairs of 
paths with average energy e in the presence of self-energy- 
only interactions. From Eq. ©, it is seen to satisfy the 
familiar Dyson equation |5J), with solution 



/^sclf 



(C) = 



1 



E„ 



(10) 



where the effective Cooperon self-energy is given by an 
integral over all momentum and frequency transfers to 
the environment, 



^sclf 



(c) 



^j{du^){dq)^ 



g,bare(^'^) 



(11) 



.,1 ■^£:,sclf 

With Zj 



q.o,barc('^'^) given in Eqs. l)A5b|l to (|A5d|l . 
Now, tlie standard way to extract the decoherence rate 
from Eq. (|10|l is to expand the self-energy in powers of 
E„ — iLo: 



vself 



(CO) 



(/:,self ^ 



_ c^,sclf ^ 



(12a) 
(12b) 



The leading "Cooperon mass" term can be identified with 
the decoherence rate, because Eq. (|9b|) yields (after per- 
forming the J (duj) integral by contour integration) 



(dq) e 



1 



(13) 



Since the {dq) integral is dominated by small g, let us 
replace g by in 7^;^^°'^ and S^^^'^, so that they can be 
pulled out of the integral. This yields 

C'"(0,i)e-^="*(*)(l 



C',^'="(0, t) 



, (14a) 
(14b) 



in which Cf^^{Q, t) is expressed in a form reminiscent of 
Eq. a free Cooperon, times the exponential of a de- 
cay function, times a factor 1 -I- S^^q'^ that renormalizes 
the overall amplitude of the Cooperon {i.e. it corresponds 
to "wave-function" renormalization, in analogy to the oc- 
curence of a finite quasiparticle weight Z in a Fermi liquid 
due to the short-time decay that is not resolved further 
by this approximation). 



Since we have to set Eq = iuj in Eq. (|12b|l and g — in 
Eq. 114b|l , it is natural to split the self energy of Eq. 11111 
into two parts, T.f^{uj) = S|';^f''i°'^(w) + Y.f^-'^{u}), chosen 
such that lyf^q'^ {lu) vanishes when Eq = ioj and q = 0. 



[This requirement is in fact fulfilled by (and was the mo- 
tivation for) the separation of Eq. ljA5b|l into two terms, 
labeled "dec" and "Z".] Thus, 7^0''" depends only on 
^sc^f,doc(^). ^gjj^g (jX5^ in Eqs. ^ and it can 

be written as follows for not too large magnetic fieldsii 
(7h/T « 1): 



09, self 

%h = 



1 f _ _ 2E° 
J (d^)(dq) 



Here the effective propagator (VV)"^^^ arising in Eq. ((T5|) 
turns out to be precisely the Pauli-principle-modified 
propagator of Eq. fI I66b|l which we conjectured by heuris- 
tic arguments in Paper I, Section IV Dl 



Ti{yy)Z - im£7(^){coth 



2T 



(16) 



-i tanh 



e — u) 



2T J 



■i tanh 



'e + u) 



2T 



The coth + tanh combination occuring in {VVyZ hm- 

its the frequency integral to \Cj\ < T, as anticipated by 
the Golden Rule discussion in Section IV CI of paper I [cf. 
Eq. {l^]. 

After performing the {dq) integral in Eq. 1)15(1 . the re- 
maining {duj) integral turns out to have an infrared di- 
vergence for quasi 1- or 2-dimensional samples. To be 
explicit, if we regularize it by hand by inserting a step- 
like cutoff function 9{\uj\ — (Dq), we obtain, for the quasi- 
d-dimensional case 



VP, self 



with pi 



duj 



2 ./i-o ^ 



tanh 



l-d/2 \ 

e— Co 



coth 



UJ 



2T 



(17a) 



2T 



— i tanh 



2T 



where 71 ~ D{e^ /haiY 



l/(27r52) and 



92 



l/(V2^7r2), 
and 73 — 



D{e^ /h(j^)~^ . For d — 2>, the integral is well-behaved 
in the limit — > 0, but not for d = 1,2. For example, 
in the quasi 1-dimensional case, the integral evaluates to 



a?,self 

7o,o 



2T 


"27i' 


1/2 


1 - 






i/2y 




. ^0 . 








V T \ 





(17b) 



which diverges for (Do ^ 0. This infrared divergence 
arises because in the present approach we have neglected 
vertex terms, which in general ensure that frequency 
transfers smaller than the inverse propagation time 1/t 
do not contribute [cf. Paper I, Section ITTi] . Thus, we 
should choose the infrared cutoff at wq — 1/i (as noted 
in Ref. ^), obtaining a time- dependent^ decay rate, 

2y27^^l/V(T^,r)^/^ wheroii 



c^.sclf 



1 



2/3 



2/3 



(18) 
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is the decoherence rate first derived by AAK^. 7^9''^^ 
grows with time, because with increasing time, the 
Cooperon becomes sensitive to more and more modes 
of the interaction propagator with increasingly smaller 
frequencies, whose contribution in Eq. (|17a|l scales like 

Alternatively, instead of cDq = the choice ujq = 
7^Q°'^ is often made, since in weak localization the time 
duration of relevant trajectories is set by the inverse deco- 
herence rate. Then Eq. I|17bl) is solved selfconsistentlyS*^, 
yielding 7^^f' ^ = (2V2/7r)2/3/r**-, with r^^ again 
given by Eq. lfTH|) . 

The decay functions for d = 1 corresponding to the 
above two choices of ujq in Eq. Ijl7b|l are, respectively 
[from Eq. (|14b|) ]: 



(2V2/7r) (i/r^r)'^' ' 



(19) 



Evidently, both equations describe decay on the same 
time scale t^^^- The second choice does not properly 
reproduce the 3/2 power law in the exponent that we ex- 
pect from Eq. ('I I44|I for F^^^{t) (a fact strongly criticized 
by Golubev and Zaikin (GZ)^). However, the first choice 
does, up to a numerical prefactor, whose precise value 
can not be expected to come out correctly here, because 
it depends on the shape of the infrared cutoff function 
(arbitrarily chosen to be a sharp step function above). 
We thus recover the classical result for the decoherence 
rate. The reason is essentially that Pauli blocking (rep- 
resented by the tanh terms in (VV)'^'^) suppresses the 
effects of quantum fluctuations (represented by the -1-1 
in coth((D/2T) = 2n{uj) + 1) with frequencies larger than 
T, as discussed in detail in Sec. IV Dl of paper I. More- 
over, we also obtain the important result that the first 
"quantum correction" to this classical result that arises 
from self-energy terms [the O ([u)o/T]^/^) correction in 
Eq. (|17b|l ]. is smaller by a factor h/y^tT). For t ^ t^, 
this is <C 1 in the regime of weak localization [cf. discus- 
sion after Eq. (ItZI)], in agreement with the conclusions of 
Vavilov and Ambegaokar.^* 



III. BETHE-SALPETER EQUATION IN THE 
POSITION-TIME DOMAIN 



The infrared divergencies mentioned above are cured 
as soon as vertex diagrams are included. However, as 
mentioned at the end of Section IH Bl and detailed in 
Appendix lA 21 frequency-entanglement then renders the 
momentum-frequency version of the Bethe-Salpeter 
equation intractable. This suggests that we try a more 
pragmatic way of finding an approximate expression for 
the full Cooperon: inspired by the insight from Paper I 
(Sec. IHI Dp that in the case of classical noise, a rather 
accurate description of the Cooperon can be obtained 
in the position-time representation by reexponentiating 



its expansion to first order in the interaction [Eq. 
of paper H], we shall try a similar approach here: we 
transcribe the Bethe-Salpeter equation to the position- 
time domain to obtain an equation for the correspond- 
ing Cooperon C'^{ri2,ti,t2) of Eq. l(5cjl . and solve this 
equation approximately with an exponential Ansatz; this 
Ansatz will turn out to yield precisely the reexponentia- 
tion of C'^'--^\ri2;ti,2t), the first-order expansion of the 
full Cooperon, in full analogy to Eq. (0. 



A. Transcription to time domain, exponential 
Ansatz 



Let us now consider the Bethe-Salpeter Eq. © for 
Cq ^^^(r^i, 1^2), i.e. with £ = s — ^^2, as needed in 



Eq. (|3b(l when ujo — 0. This equation can be transcribed, 
using Eq. lf5c|) (with there replaced by e) to the form 

{-DVl^ + dt, + 7h) C'iri2;h,t2) = S{ri2) <5(ii2) + 

dr^dtidt'^ Sf;^i;'(ri4;ti,t4)C"'(r42;i44',i24') , (20) 



where the self-energy in the energy/position/times rep- 
resentation is defined by: 

Sf;if(ri4;ti,t4) = j {dq){d^l^){d^2){.dni) (21) 
^ ^[qr,,~n^t,+n,u-t'^(n^~n,)]^e-\ii2^^^^^^~^ 

Before trying to solve Eq. H2U|I . let us get a feeling for 
the structure of this equation, by calculating the zeroth 
and first order terms of C^(ri2; ti, ^2) in an expansion in 
powers of the interaction propagator [i.e. Share)- To this 
end, we use the fact that 

{-DVl + dt+ 1h) C\r, t) - 5{r) 5{t) , (22) 

iterate Eq. (|20|) once, and replace Sfuii by Share [given 
by Eqs. IjASp ] on its right-hand side. We find 

C'^{ri2-MM) = C\ri2M2) + C^^'>'{r^2-MM) + ■ ■ ■ , 

(23a) 

where C''-^''^(ri2; ti, t2) has just the structure discussed in 
Paper I, Sec. IHI CI describing propagation from (r2, t2) to 
(ri,ti), with interaction vertices along the way at points 
(r4,t4) and {r^.tz): 

&^^'{ri2;h,t2) = J dr3dt3dr^dtiC°iri3,ti3) (23b) 

X^baro(^34;i3,i4)C'°(r42,t42) , 
SLre('^34;t3,i4) = / < S^f, (r34 ; ^3 , ^4) (23c) 



= J (d(7)(df73)((if74)e*['"'^^""^*^+"***l 

xs^l'^^(^^3,^^4). (23d) 
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Let us now construct an approximate solution of the 
Bethe-Salpeter equation l|2U|) . by making an exponential 
Ansatz of the following form: 

C%n2;h,h) = C'"(ri2,ti2)e-^^(''-^*^'*=) . (24) 
The decay function F"^ is needed only for ti2 > (since 



C'^{ri2,ti2) vanishes otherwise), and is required to obey 
the initial condition ^^(0, t2, t2) = for all t2- The 
Ansatz (I24II solves Eq. ()2nil exactly, provided that the 
decay function F'^ satisfies the equation 



C°(ri2,ii2) 



dt, - DVl^ - 2D 



VnC"(ri2,^i2)- 
C'°(ri2,ii2) 



{ri2 ;ti,t2)-D [Vr, (ri2 ; ti , i2) 



dr, dU dt'^ (ri4; t,,U) C\r,2M2) 



(25) 



B. Evaluation of the decay function F^{t) 

Let us now evaluate the decay function F'^ (t) explicitly; 
after three simplifying approximations, we shall find that 
it reproduces the function F^^^^{t) of Eq. fl l^ . 

Our first simplifying approximation is as follows: in- 
stead of trying to solve Eq. (|25|l in general, we shall be 
content to determine the decay function F^ only to linear 
order in the self-energy, in accord with the fact that we 
"linearize" the latter by replacing the full self energy by 
the bare one. (Including nonlinear contributions would 
add terms that are smaller than those kept by powers 
of the small parameter 1/g ) To this end, it suffices 
to hnearize Eq. (|2SJ) in F^ , by dropping the {VnF'^Y 
term on the left-hand side and the exponential factor 
e~^^i2-^i2\ on the right-hand side, and replacing Efun 
by Sbaro- One readily finds that the resulting linearized 
equation is solved by 



F^ri2;ti,t2 



C0(ri2,ti2) 



(26) 



where C^^^^ is given by Eq. (|23b(l . Thus, the expansion 
of C" [Eq. (|23Il] to first order in F" reproduces Eq. (|23a|l . 
as it should, and conversely, Eq. turns out to be 
nothing but the reexponentiated version of Eq. (|28a|l . 
Our explicit solution of the Bethe-Salpeter equation, to 
linear-in-S accuracy in the exponent, thus very nicely 
confirms the heuristic analysis presented in Section IlII Dl 
of Paper I in favor of reexponentiation strategies. 

The second approximation is necessitated by the first: 
upon comparing with the structure of the self-energy- 
only solution [Eq. H14a|) ]. and following the discussion 
before Eq. (|15|l . we recognize that effectively only a 
part of C^^^*^ may be reexponentiated (note that this 
remark would be irrelevant if we were able to find the 
exact F'^). Therefore, when evaluating C^-^"^^ explic- 
itly from Eq. ^B, we insert E^are = ^HL + ^IZl 
[Eqs. (jASap ] into Eq. (|23d|l . but for the self-energy term 



^barc l|A5b|) ] we retain only the "decoherence" con- 

Eq. |J5d)l ] 



"'bare 

tribution S 



clt.dc 
bare 



[Eq. (|A5c|l ]. because Y.' 



■self,Z 
bare 



contributes only to the renormalization of the overall 
amplitude of the Cooperon (and in any case its time- 
dependence for long times turns out to be weaker than 
that arising from Sf^'^!^*"^, as is checked explicitly in Ap- 



pendix [^21) ■ In other words, we write 



^soif Z' ^^'^ drop the second term. The resulting expres- 



dec 



sion for C, 



C^'J:{ri2;h,t2) - ^ / {dq){duj)idq){dLu)e^'^^^^ e"^"*- 



reads 



X C°(c.)E«^(..,^)C^(u;) 



(27) 



+ e" 



*12 nOf, , -^ y^e.vert 



(c^,cS)C°(a; + c^) 



[The quickest way to arrive at Eq. H27|l is from the sec- 
ond term of Eq. with (fig, ^^2) ^ 27r(5(r223)C^(r22) 
on the right-hand side, and Sfun — > Sbare, given by 
Eqs. 1IX5I1 .] 

Our third approximation for evaluating the first or- 
der decoherence correction to the Cooperon (and thus 
the decay function) consists in retaining only its domi- 
nating long-time behavior, for Tti2 S> 1. In this limit, 
terms of order ui/T are <C 1 and may be neglected (they 
produce subleading contributions for F'^{t), as is checked 
explicitly in App. IC This allows us to keep only the 

= component of the effective environmental propa- 
gator Cf^^^{Lu) [Eq. (|A5f|l ] . which is contained in both 

^7ZT and see Eqs. and More 

formally, after substituting the latter two equations for 
the S's occuring in Eq. (|27|l and symmetrizing the in- 
tegrand w.r.t. UJ <-> —Lu, we Taylor-expand £^^'^^(0;) in 
powers of lo and represent uj as «9ti2 under the Fourier 



7 



integral, thereby bringing Eq. (|27|l into the form 

n <J 



ri=0 



xC'^(ri2,ti2)7'(;™,,^,,,)(g-c^). (28) 



Its ingredients are defined as follows: 



2n\ 



(29) 



— crw . _ _ . 

>^(r-i2;ti,t2)(9,t^) = 



2 {d^){dq)^ 



rti rt3 

2 dts dt^ 

Jt2 Jt2 



C^(c^) Cy,{uj - Lu) C'^iu;) - e^^*- C°{u; - u,) C^_,-(l.) C^(a; + u,) (30a) 



-P(riltl2) (9^*34) 
I 



(30b) 



Eq. (j30b|l follows by transforming the free Cooperons in 
Eq. H30a(l to the time domain, and recognizing the (dq) 
integral of the resulting expression to contain the object 



-P(r™,tl2)('?'^34) - 



C0(ri2,ti2) 



*34 ' 
*12, 



(31) 



This quantity is the Fourier transform of the probabil- 
ity density P(o,t) (''': ^0 foi' intermediate portion of a 
random walk to cover the distance r' in the time t', un- 
der the condition that the total walk is closed, returning 
to the starting point r = after a total time t. It was 
introduced in paper I [Eqs. fl ES|) and (I|25J)] as central 
ingredient for averaging the effective action of the influ- 
ence functional derived there over pairs of time-reversed, 
closed random walks. 

We are now in a position to write down an explicit 
expression for the decay function -F^(ri2; ti, t2)- Writing 
= Er=o ^{n}' find from Eqs. l|2HIl and Eq. 1^)) : 



1 



C0(ri2,ti2) 



n (32) 



x/(dg-)C"^(ri2,ti2)P(?,:.*,,)(g-i34) 



We henceforth set ri2 = and ii ~ — ^2 = ^t, as 
required for calculating the conductivity [cf. Eq. (jSbl) ]. 
and write F^^-^it) = F^^^{0; ^t, ~^t). We shall dis- 
cuss only the leading term ^(q), since the -F'(n>o) terms 
give subleading contributions. (This is illustrated in 
Appendix IC II for i^(i).) For n = 0, the correlator 
needed in Eq. 122) reduces [via Eqs. if^ and HA5f]l ] to 
the Pauli-principle-modified noise correlator of Eq. 
^efn) qi'^) ~ After symmetrizing the range of 



t /2 

the i4 integral to be J_^/2'^^'t^ ^^'^ setting ti —t^ in 



the e'"*34 term of Eq. we obtain 

" 2 

^(07)(9" l^34|)-P('o7)(9" 1^341) 



dh j dU I (dq) / (dw) e"*"*-'^ (33) 
pp 



This result is identical to the function i^P.P,(t) whose form 
was conjectured by heuristic arguments in Section fVDI of 
paper I, namely Eq. fI I65|) [with SP therein given by the 
"closed random walk" version of Eq. fI I32al) ]. Thus, we 
reach the main conclusion of paper II: the heuristic way 
of introducing Pauli blocking into an influence functional 
approach in Paper I, Section IV Dl (and, by implication, 
also the more formal analysis oTFefB My cons^stent 
with the present diagrammatic Bethe-Salpeter approach. 

To calculate the energy-averaged version of the 
Cooperon, ((7^(0, i))^, as needed in Eq. l|5all . we need 

the energy-average of e"^'"'*-*'. This was done in great 
detail in Paper I, Sec. I VI K\ for d ~ 1,2,3, so we shall 
quote only the result for d = \ here: 



(34) 



(i?fo)W>e = F^Ut)h + 0[{tT)-'/']] . (35) 



Here F.^Ut) = iV^/m/ 



AAK^3/2 

' <p,l I I 



with 



given by 



-t/2 



Eq. ©, is the result for the decay function F^l^{t) ob- 
tained in paper I [Eq. fI I44(l ] for classical white Nyquist 
noise. Eq. H35|l states that the leading quantum correc- 
tion to the F^l^{t) is of order (Tt)^^/^ oc 5^^^^, i.e. small 
in the regime where weak localization theory is applica- 
ble. [The numerical prefactor of this term was evaluated 
exphcitly in paper I, see Eq. (|73a|) ]. 



C. Comparison with magnetoconductivity of A AG 

As a final check of our Bethe-Salpeter analysis, let us 
use it to directly calculate 5a^^'''^\ the first term in an 
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expansion of the weak localization conductivity in pow- 
ers of the interaction. Inserting C^l!^' {0; ^t, —^t) from 
Eq. (EH) into Eqs. ^ and Eq. ^ to obtain Ja^^^^), 
the / dt integral produces a S{lu) that sets a; = 0; af- 
ter some obvious substitutions [from Eqs. l|A5cp . l)A5e|l 
and we readily find that the resulting expression 

for (5(7^'"'^' is given precisely by Eq. fI I76() . which, as 
mentioned previously, agrees with Eq. (4.5) of AAG^. 



D. Plausibility arguments for exponential Ansatz 

To end this section, some remarks on the adequacy 
of our exponential Ansatz are in order. First, if an ex- 
act solution for Eq. for F^(ri2; ^i, could be found, 
the exponential Ansatz H24|l for C^{ri2;ti, would yield 
the exact expression for the Cooperon. Of course, how- 
ever, it was necessary to make approximations in solv- 
ing Eq. 1)25(1 . and once these have been made, one might 
question whether the exponential Ansatz adequately cap- 
tures the important physics. For example, one might 
consider functional forms of the type Ad{t)e''^'^^^\ as 
discussed by Golubev and Zaikin2i, where the prefactor 
Ad{t) has a non-trivial time dependence different from 
C". [With such an Ansatz, it would not be possible to 
determine Ad{t) and Fd{t) from a first-order calculation 
of the Cooperon, since it would be unclear how to sep- 
arate the contributions of Ad{t) and Fd{t) to C*^ in or- 
der to decide which part has to be reexponentiated and 
which part should stay in the prefactor.] Indeed, GZ have 
arguedSithat the final expression for Fd{t) after averaging 
over diffusive paths contains only coth[cZ>/2r] factors and 
no tanh[(e ^ tD)/2T] factors, and that the latter instead 
only contribute to the prefactor Ad{t), in such a way that 
an expansion of Ad{t)e~^'^'^^^ to first order in the inter- 
action propagator correctly reproduces the combinations 
coth-|-tanh occuring in C^. In our language, that would 
correspond to reexponentiating only contributions from 
S^, while attributing all contributions from to the 
prefactor Ad{t). 

However, there are several strong arguments against 
such a procedure. Firstly, the diagrams for and al- 
ways occur in matching pairs, generating a series of prod- 
nets of the type C^") ([S^ + S«]C(o) . . . C^o) ([E^ -h S^]C(o). 
An Ansatz assigning to the prefactor only and S/ 
to the exponent would disrespect the structure of this 
series. Secondly, the structure of the Ansatz should be 
sufficiently general that it holds for any dimensionality, 
d = 1,2,3. But for d = 3, the complicated Bethe- 
Salpeter analysis is not necessary and the Dyson equation 
sufficient, because the self-energy contributions 
[Eq. 1)11(1 ] are infrared-convergent by themselves, without 
the need for vertex corrections; in this case, the Cooperon 

decay is indeed purely exponential, e"'''" " *, where the 
decoherence rate is given by 7^_f " of Eq. (|T73, which 
evidently does contain the combination coth((D/2r) -|- 



tanh(£ — a;)/2T. The fact that this combination shows up 
in [t) implies that it should also show up in 2 and 1 di- 
mensions for the decay functions i^i,2(i)- Thirdly, in the 
limit e ^ T, it is general consensus that the decoherence 
rate is simply given by the inelastic rate, ~ e'^^-^t, 

and indeed this result was recovered from our theory in 
paper I [Sec. lVI C] : but this is possible only if F'^{t) con- 
tains tanh[(e — uj/2T)] functions, since they are the only 
way in which the energy dependence enters the theory. 
And finally, the fact that the combination coth + tanh 
occurs in the effective action of the influence functional, 
i.e. in the exponent, was derived by GZ themselves^ us- 
ing their influence functional approach (the fact that the 
tanh-terms dropped out of their flnal results for the de- 
coherence rate is only due to their neglect of recoil, as 
shown in Ref. 0). 

IV. CONCLUSIONS 

In papers I and II, we have shown how the combined 
effects of quantum noise and the Pauli principle can be 
incorporated into a calculation of decoherence rate of in- 
teracting electrons in disordered metals. To this end, we 
used both an influence functional formulation and stan- 
dard diagrammatic methods, obtaining identical results 
with both methods. The influence functional approach is 
perhaps more intuitively transparent: it is formulated in 
the position-time domain, where we have intuition about 
the behavior of diffusive trajectories, and shows very 
nicely how for quantum noise the contribution to deco- 
herence that arises from spontaneous emission gets can- 
celled by Pauli blocking at sufficiently low temperatures. 
Thus, we find that as long as the condition Tr,^ ^ 1 
holds (which characterizes the regime of weak localiza- 
tion), the decoherence rate decreases without saturation 
as the temperature is lowered towards zero. The fact that 
Golubev and Zaikin obtain a saturation was identified to 
be due to their neglect of recoil. 

However, paper I does rely on heuristic arguments in 
the way Pauli blocking is introduced. Corroborating the 
correctness of these arguments was the purpose of paper 
II; and indeed, setting up a Bethe-Salpeter and solving it 
by an exponential Ansatz, we recovered the decay func- 
tion found in paper I. Moreover, we identified several 
correction terms that do not arise in the influence func- 
tional approach {F^^\ ^sc\iz) ^^'^ showed them to be 
neghgble. 

Apart from clarifying the fundamentally important in- 
terplay between spontaneous emission and Pauli block- 
ing, our calculation of the decoherence rate has the merit 
of being free from any infrared or ultraviolet divergen- 
cies: in the leading terms that govern the decoherence 
rate, all necessary cutoffs arise naturally from within 
our formalism, and do not need to be inserted by hand 
(whereas AAK did need to insert an UV cutoff by hand 
for d = 2,3). This has enabled us to obtain a number of 
new results. Firstly, our more accurate treatment of the 
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regime of large frequency transfers {Co ~ T) has allowed 
us to calculate explicitly the leading quantum corrections 
to the results of AKK for T^p [Eqs. fI I75(l ]. finding them to 
be small in 1/Tt^. Secondly, by explicitely keeping track 
of the energy dependence of the propagation energy e of 
the diffusing electrons, we were also able to discuss in de- 
tail the energy-dependence of the decoherence rate, also 
for energies higher than the temperature Eqs. (I I85|I . 
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Appendix A: DIAGRAMMATIC DERIVATION OF BETHE-SALPETER EQUATION 



This appendix diagrammatically specifies the Bethe- 
Salpeter equation governing the Cooperon, and gives ex- 
plicit expressions for the full and bare Cooperon self- 
energies occuring therein. 

1. Full Bethe-Salpeter equation 

The diagrams specifying the Bethe-Salpeter equation 
for the full Cooperon, first written down in Ref. i3, are 
shown in Fig. |2| The feature that distinguishes "weak- 
localization" from "interaction" corrections to the con- 
ductivity is that for the former each current vertex is 
attached to a retarded and an advanced electron propa- 
gator, G^jG^, whereas for the latter, one or both cur- 



rent vertices are connected to two retarded or two ad- 
vanced electron lines, G^jG^ or G^jG^. Hence, in set- 
ting up the Bethe-Salpeter equation, only those types 
of diagrams have been included for which the upper (or 
lower) lines entering and leaving the Cooperon are both 
retarded (or advanced) electron propagators. By adopt- 
ing a pure ladder structure in Fig. |2Ia), diagrams con- 
taining overlapping interaction lines have been dropped, 
but these are known^^ to be smaller than those included 
by at least a factor of (fep^oi)^^ ^ 1 (where Zoi is the 
mean free path). 

Fig. |2Ia) translates into the following equation: 



1 ^ 



—c^jni,n2) = 27r^(f7i2) + 



Tel 



{Eq- ifll) 

Tel 



c'q{n^,n2)+ /(dr!3)Eff,ii(f)i,f]3)c^(f^3,f^2). (Ai) 



This equation is equivalent to Eq. JHl of the main text. 

The "Cooperon self-energy" j^jj occuring herein is defined by 



v^f self o -\ I ■sT'f ,/+-R,vcrt o -\ 

Sq.g-fuU (f^l'^3,w) +Sg:,-,fun (f^l,f^3,w) 



(A2a) 



where the subscript S + indicates a sum S -I- S , and the self-energy and vertex contributions to S are given 
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Figure 2; (a) Diagrammatic depiction of the Bethe-Salpeter equation Eq. ijAl^ for the Cooperon Cq{Q.\,Q. 



The contributions to the full self-energy S f^ll (fii.f^a) of Eq. Eli, namely (b): E^'^^/^if "(J^i 



The superscripts I /R on E indicate which type of interaction 



(which contributes to Ej, 



■£,//«, vort 



only); and (d) Elj^^f^n 



(b), (c) and (d): 
"9.?,] 



propagator occurs: E (generated by Si in GZ's approach-^) contains C , together with Q and (J on the upper and lower 
contours, respectively; E (generated in GZ's approach by Sr) contains a sum of two types of terms, featuring the combinations 
together with and on the upper and lower contours, or together with and on the upper and lower contours, 



■//-R.vert 
full 



occur both without and with bracketed labels for the interaction and electron 



respectively. The diagrams in (d) for E 
propagators, e.g. R or [A) on C, to distinguish contributions in which the electron Keldysh Green's function occurs on the 
upper or lower electron line, respectively. The vertex to which the electron Keldysh Green's function is attached is always 
indicated by a double dot (adopting a convention used in-), (e) Diagrammatic conventions; the dressed interaction vertex is 
denoted by Vq{Lu) — \V'^{ijj)]/Tc\, where is the bare diffuson propagator. 
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diagrammatically by Figs.jSfb) to|2fd). These lead to the following expressions (here nf = £±^njj 



-£,_R,solf 



S,:,-fun (f^i>f^3,cS) = tanh[(f7+ - lu)/2T] ^zC^iCu) C^J"" (n^ - ^,^3 - cS) 



- tanh[(f77 - c:;)/2T] ii£^(cj) C^_|"'(f7i + cj, f^s + cD) 



(A2b) 
(A2c) 



^l',j!'Hikami(^l'^) 



^ tanh[(f7+ - Cu)/2T] ^^Cf{Cu) [V'^iCo)]' ([C°(f^i)]-' + [V'^iCo)]- 

-tanh[(f7r - u^)/2T] (p^(^^l)]"' + [^°(-^)]~' 



(A2d) 



z/:f(c^) 



1 

"2 "'-'9 



(r!i + CD, rig - + cr'fisix + cs, rig - c:-) 



(A2e) 



[tanh[r!+/2r] £f (cD) - tanhi^ig /2r] C\Xf (Six + 1^3 - 



(A2f) 



+ [tanh[(r!+ - cD)/2T] - tanh[(r!^ - w)/2r] £^(w)] c\_f{^x + cj, fia - (i) 



The terms in Eq. l|A2dp . with their characteristic depen- 
dence on (pO)2[(c°)"Vf£")"^], stem from the Hikami- 
box contributions of Fig. Ejc). The ingredients entering 
the above equations are given by 



'CT(cD) 



4 



2i coth(<I./2T)Im[/:|^(a;)] 

mm* 



q 

iu) 



q 

2vE°^ 



(A3a) 
(A3b) 

(A3c) 
(A3d) 



In Eq. (|A3b|) we have taken the usual "unitary limit", 
which is relevant in the limit of small frequencies and 
momenta; the more general expression isi& 



2uEl + (RO ~ iCo)/V^''^ ' 



(A4) 



where 



(d) 



Id' 



r e 



-iqr( 2 



(e /r) is the Fourier trans- 
form of the Coulomb potential in d effective dimen- 



sions: ' — 
a?e^ \n{(fa?). 



^q 



(2) _ 



ae2 27r/|g|, and 



(1) _ 



2. Bare Self-Energies 

Since the self-energies Efuii of Eqs. (|A2(I are propor- 
tional to the Cooperon C, the Bethe-Salpeter equation 
is non- linear in C. Solving it in its full glory thus 
seems hardly feasible. Thus, we shall "linearize" it by 
making the replacement C^{^i,Vtz) — > 27r(5(ni3) C°(ni) 
for every occurence of the full Cooperon in the self-energy 
terms. The resulting bare self-energies can be written in 
the form 



- / {dq) / {duj) 2TT 



(A5a) 
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where r2i4 = fJi 



(j,g,barc V 



E£^,sclf,dcc/ - \ 



and rii4 = ^(rii + ^4), with the following ingredients: 



(J, bare q,g,barc 



(A5b) 
(A5c) 



E£,sc\f,Z / - \ 
9,g-barc('^'^) 



tanh[(£ - c^)/2T] izZf (cS) ' ( [C°(.;)] + [pO(c^)] - C0_,(^ + u) 



(A5d) 



- tanh[(£: -u- u;)/2T] ^tCfi^) [V^i-co)] ' ([C°(c.)] ' + [V'.i^cu)] ' ) C'^_^{cu - Cu) 



E 5, vert / -\ 



"'g,q,barc 



2C°_,(c^)£t%-(<:^), 



= coth[(I>/2r] Im[/:f ((!>)] - tanh[(f ~ u})/2T] \iC^{Ld) + tanh[(f - uj - uj)/2T] \iC^{Cj) 



(A5e) 
(A5f) 



We have split the self-energy contribution [stemming 
from Eqs. l|X2b|) plus l|J2c|) ] into two terms, S^"";'^'"' + 
^bare^' chosen such that the Hikami-box contributions 
are fully contained in '^^^^^'^ , and that both 2^'!^^^^'^°'^ 
and SJJ^^j.* are proportional to the same combination of 
propagators, C^'^'^ [Eq. ljA5t]l ] . a feature that consid- 
erably simplifies the analysis in the main text. To 
achieve this, the terms in Eq. (|A5dl) that are propor- 
tional to tanh • Cg_ g{uj zk iD) were added in Eq. ljA5c|) 
and subtracted in (|A5dp . respectively. This addition- 
subtraction trick amounts to " replacing" the Hikami-box 
contribution toSf°'^^'^'^'^ by "replacement terms" [those 
added to Eq. (|A5c(l ] that (i) have a simpler, more con- 
venient structure (since proportional to C° instead of 

(P0)^[(C°)"^ + {V°y^]), but (ii) nevertheless have the 
same leading infrared and ultraviolet behavior, in the 
sense that the difference between the Hikami-box and the 
replacement terms, namely E'^^'f'^^ generates only sub- 
leading contributions to the long-time behavior of the 
Cooperon (as explained below). The leading contribu- 
tion comes from S^^oif.dcc g^^^ S™''*, because both are 
proportional to the effective propagator Cf[^^^{Oj), whose 
coth[(I'/2T] term at small frequencies lj makes the domi- 
nant contribution [in contrast, S'^'^'f'^ lacks such a term]. 

Although the dominant contribution comes from low 
frequencies, 2'^°'^''^°'^ -I- E^'^''' contain no infrared diver- 
gence, since for w — > 0, their contributions cancel each 
other, as is clear directly from Eqs. (|A5cp and l|A5e|) (or 
by inserting them into Eqs. (|A6b|) and (jA6cp below, or 
Eq. (|27|l of the main text). Moroever, the dominant con- 
tribution from 5]^°'^''^°'^ + S™''* also contains no ultravio- 
let divergencies, since the effective propagator £f^^^{cj) 
evidently vanishes exponentially in the limit uj 3> S^T, 
hence the / (dOj) integrals over both 2]^'^''^'^|^°'^(a;, w) and 

^g'lj baro('^' Separately free from ultraviolet diver- 
gencies. 

Concerning the contribution from Y^^'^^^'^ ^ there are sev- 
eral ways to convince oneself that its contribution to the 
Cooperon decay function F(t) is subleading. Firstly, an 



explicit calculation, performed in Appendix IC 21 shows 
that its contribution to C'^^-''^, the first-order expansion 
of the Cooperon in powers of the interaction propaga- 
tor, depends much more weakly on propagation time t 
than that from 5]''°'^''^°'^ -|- S™''*; for example, for quasi 1- 
dimensions it scales with i^^/^ [cf. Eq. (|C15p ]. compared 
to the t^ of the leading terms (^"(O, t)F^{t). 

An alternative, and more simple, argument goes as fol- 
lows: according to the short-cut "self-energy-diagrams- 
only" approach of Section III CI the decoherence rate is 
given by7,^b ofEq. ^ [see also Ijl4b|l ] , for which we 
have to set Eq = icu and take g — > 0. Now, in this limit 
'^Tq '^ i^) 'Vanishes identically, regardless of the form of 
the interaction propagator Cq ifu), and hence does not 
contribute to 7^0°'^ at all. Actually, an even stronger 
statement can be made if the general form (|A4|) for the 
interaction propagator is specialized to the the so-called 
"unitary limit" of Eq. (|A3b(l (as is usually done anyway). 
In that case, both the Hikami-box terms and the "re- 
placement terms" separately vanish in the limit Eq = icu 
and 9 — > 0. Thus, for the "unitary limit" of the in- 
teraction propagator, Hikami-box contributions actually 
do not contribute to the decoherence rate at allS^. It is 
for this (somewhat fortuitous) reason that the influence 
functional theory of decoherence developed in Ref. was 
able to correctly obtain the Cooperon decay function, de- 
spite of the fact that it did not include any Hikami-box 
contributions. 



3. Expansion to 2nd order in the interaction 

The object needed on the right-hand side of Eq. lf5c|) 

for C^(ri2, ii, t2) is ^ ^ (fii, f22), which is determined 
by the Bethe-Salpeter equation Eq. with £ = e — ^il2 
(with loq — here). It is instructive to consider the first 
few terms that are obtained upon iterating this equa- 
tion, while using the bare self-energies of Eq. (|A5a|) . (Of 
course, as soon as we go beyond first order, we should not 
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use the bare self-energy, but the full one, which should 
be calculated iteratively order by order, too; we shall re- 
frain from doing so here, since our intention is merely to 
illustrate the general structure of the terms arising in sec- 

I 

C(o) 



These expressions are useful for illustrating two impor- 
tant general points. Firstly, the expansion l|A6|) allows 
us to confirm explicitly a fact well-known to practition- 
ers of diagrammatic perturbation theory, namely that the 
calculation of the Cooperon is free of ultraviolet divergen- 
cies. This fact was implicitly challenged by GZ, whose 
conclusion of a finite decoherence rate at zero temper- 
ature stems from the occurence of an ultraviolet diver- 
gence in their expression for the decoherence rate. (GZ's 
expression for 7^^ for d — I [Eq. (76) for Ref. ^] has 
the form of our Eq. H17a(l . but without the tanh-term, 
whence they introduced an upper cut (Dmax — I/toI in 
the frequency integral there; the relation of their work 
to ours is discussed in more detail in Paper I, Sec. IVIII ) 
However, it is straightforward to check [using Eqs. IjASfl ] 
that the perturbative expansion l)A6|l generates no ultra- 
violet divergencies when used to calculate that version of 
the Cooperon governing the conductivity, namely C^^'l^ 
of Eq. H3bp : the reason is simply that both S'^"'^-'^®'^ and 
S™'* are proportional to the propagator 'C^^'^^(a)), which 
serves as ultraviolet cutoff at a) ~ T. [The contribution 
from is subdominant, as mentioned above; in fact, 

'^q 'q^hiirci^ ' ^) ^ in the limit of large oj, and its lead- 
ing nonzero contribution turns out to be UV-convergent 
if the general expression ljA4|l is used for the interaction 
propagator, instead of its small- frequency approximation 
iTOhll .] 



ond and higher orders, not to evaluate them explicitly.) 

Cg~^^^ni,n2) = C(°) -KC(i) -HC(2), and evaluating the 
result for fli^2 = ti-" T t^, we obtain: 



(A6a) 
(A6b) 



(A6c) 



Secondly, we note that the frequency arguments uJi get 
more and more "entangled" from order to order in pertur- 
bation theory, i.e., they occur in increasingly complicated 
combinations as arguments of 5]seif/vert^ because the ver- 
tex diagrams cause a proliferation of frequency transfers 
between the upper and lower Cooperon lines. In n-th 
order, the generic structure will be 

C(")~n /(dcS.)C(°^(-)S(-,c^i)...S(,(ii)C(")(.), (A7) 

where (•) stands for combinations of frequency arguments 
that can contain any number of u)j 's. Due to this entan- 
glement, a direct solution of the Bethe-Salpeter equation 
(jHJ is intractable, and further approximations are needed, 
that somehow "factorize" the entangled frequency inte- 
grals and thereby truncate the proliferation of frequency 
transfers. 

A natural truncation scheme would be to retain the 
frequency transfer ujj generated by a given vertex line 
only in the corresponding vertex function I]^'^"*(ti;, (Dj), 
and to neglect it everywhere else in the diagram. As 
a result, it would again become possible to associate a 
definite frequency label with the upper and lower elec- 
tron lines of the Cooperon (say e and e — u). In fact, 
such an approximation was in effect adopted in the in- 
tegral functional approach of Ref. and paper I (which 
both implicitly also took the "long-time limit" a; = 0, 



^ C°(a;)5(22), 

= i|(dg-0(d<ii)|.5(2i)C°(a;)S^;|f,„Jc^,cDi)C^ 

+ 5(25 - 2cSi)C°(c. - cDi) S-™;;*,,,(c.,cSi)C°(u; + cDi)| , 
= ^/ {dq,){duo,){dq2){du2)\^5{2:o)Cl^^^ 

+ S{2^ ~ 2u,) C°(c^ - c^i) S-_;-;:^_(^, c^i) C°(c. + Co,) S5;|f,_(^ + Co.^u^^) C°(a; + Co,) 
+ 5{2do - 2C02) C°(c. - c^2) - c^2, c^i) C°(a; - u^2)^'^:Lr.i^ , ^2) C^(^ + ^2) 

+ 5{2ui~ 2Cdi - 2012) C\{lo -Cbx- C02) 5^^^"bare*(^ ~ 0)2, <^i) C°(cj + Cdi- C02) 

>^'^q2%roi^ -h Wi, W2) C°(W Wl t;)2) > 
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for reasons explained in the last paragraph of Sec. I VII XI 
of paper I). Such a procedure can be justified as fol- 
lows: The only reason for incorporating the (frequency- 
proliferating) vertex terms in the first place, is to cure the 
infrared divergencies arising from the self-energy terms, 
which are thereby cut off at frequencies w ~ 1/t [this 
is perhaps seen most clearly from the [l — sin(a;t)/((Dt)] 
factor in Eq. fI I41|l ]. For larger frequencies (D > 1/t, 
the contribution of vertex diagrams is always subleading 
compared to that of the matching self-energy diagrams, 
and hence can be neglected without affecting the leading 
behavior of the decay function Fd{t). Thus, it suffices to 
treat the a)j -dependence associated with the frequency 
transfer between the forward and backward contours in 
the vertex part of E(-, (Dj) explicitly only within this par- 
ticular factor [i.e. in the interaction propagator, associ- 
ated coth + tanh functions, and associated Cooperons of 
S(-,tZ'j)]. Since the accociated contribution is dominated 
by frequencies w ~ 1/t, which are small in the long- 
time limit, all other factors C(-) and 'E[-,CdijLj) of the di- 
agram to which this (Dj-dependence has propagated may 
be Taylor-expanded in tUj. Moreover, only the zeroth- 
order terms of this Taylor expansion need to be retained, 
since the others contain higher powers of Wj ~ 1/t, and 
hence produce contributions with a subleading time de- 
pendence [as illustrated in App. IC II where such an ex- 
pansion is carried out explicitly in a very similar context] . 

Having clarified that a truncation scheme is justified 
in principle in the long-time limit, we have to implement 
one in practice, in such a way that the leading terms 
are not affected. In the present context, the simplest 
version of such a truncation scheme would be to replace 
the 5(rii4 -I- 2uj) in Eq. l|A5a|l by 6{fli4), thereby render- 
ing the entire equation for bare(^i' ^4) proportional to 
(5(rii4). As a result, the arguments of the "self-energy- 
diagram-only" discussion in Section III CI would apply: 
the Bethe-Salpter equation could then be simplified to 
a Dyson-type equation, whose self-energy would be given 
by an expression analogous to Eq. Hll|) . but now includ- 
ing a vertex contribution: 



1 



EE - / {dco)idq) 



q,g,bare 



■^e,VGrt 
q,q,havc 



(cj,cj) . 



However, note that it would now not be possible to calcu- 
late the decoherence rate 7^0 according to Eq. H12bp by 
setting Eg = iuj and q = in because the factor 

Cq_g{uj) contained in cj) would then yield an 

infrared divergence for g ^ 0. 

To avoid this problem, a version of the calculation 
has to be found in which the condition Eq — ilu is 
avoided, and the variables q and to are integrated over in- 
stead. In Section UTTI of the main text this is achieved by 
transcribing the Bethc-Salpeter equation to the position- 
time domain, and solving it with an exponential Ansatz, 
Cin2;h,t2) ~ COelc-'^Vc"], where Cf^) (ria; ti, ts) is an 
appropriately Fourier-transformed version of Eq. (jA6bl) 
involving J{dq){dw) integrals, precisely as desired. This 



is a factorization approximation, in the sense that a pro- 
liferation of entangled frequencies is avoided by approx- 
imating the n-th order contribution to the Cooperon by 

Note that in this scheme, the frequency transfer be- 
tween forward and backward lines generated by the ver- 
tex terms is treated exactly in the first order terms needed 
for C^^^; the factorization approximation sets in only in 
second and higher orders. Treating the first order terms 
exactly is the best one can do in our reexponentiation-of- 
C'-^-' scheme, since in the latter, an accurate treatment of 
effects occuring only in second or higher order is beyond 
the accuracy of the method. The accumulation of energy 
transfers is such an effect, but fortunately it produces 
corrections that are only subleading in time, as argued 
above. 



Appendix B: THE BETHE-SALPETER 
EQUATION YIELDS THE EXACT COOPERON 
IN THE CLASSICAL LIMIT 

In general, solving a Bethe-Salpeter equation starting 
from a self-energy calculated only to lowest order in the 
interaction does not provide the exact solution of the 
initial problem. This remains true even when the self- 
energy is treated self-consistently (i.e. inserting the full 
propagators into the diagram for S, as we have done for 
J^fuu)- Nevertheless, in the following we shall demon- 
strate that for d — 1 and classical white Nyquist noise, 
the exact solution of the Bethe-Salpeter equation (|20|) 
fully reproduces the exact results for the Cooperon de- 
rived by AAKi, implying that our Bethe-Salpeter equa- 
tion itself is exact for this type of noise. The reason for 
this may be traced back to the special properties of the 
white Nyquist noise interaction propagator, as will be 
explained below. Thus, non-exact results obtained from 
our Bethe-Salpeter equation for d = 1 and white noise 
are entirely due to approximations involved in construct- 
ing a solution, such as the "re-exponentation" of the first 
order result. (Actually, the deviations resulting from the 
latter approximation produces are quantitatively rather 
small, as demonstrated in Paper I [Section IlII Ej .) 

We start from the Bethe-Salpeter equation, Eq. H2U|I . 
using the full (not bare!) self-energy Efuii, given by 
Eqs. H21|) and ljA2|l . The latter simplifies considerably 
for classical white Nyquist noise, described by setting 
£f/^((I>) ^ and (w) ^ Cf = T/{vDf). The 

first replacement implies that the self energy diagrams. 



and hence also the Cooperon €^{0.1,0,2), no longer de- 
pend on E, so that this argument will be dropped hence- 
forth. The second replacement results in an interaction 
propagator that is not only independent of e, oj but also 
of uj, i.e. white in frequency transfers, implying that it 
becomes a (5-function in the time domain. The resulting 
self energy has the form 

Ef„*if(ri4;ii,t4) = W'Jru)[5{tu) - 5{hi)\5{ti,), (Bl) 
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Figure 3: Schematic illustration of the time-evolution in 
the differential equation for the Cooperon: AAK's scheme 
(top) evolves both endpoints symmetrically, while the full 
Bethe-Salpeter equation (bottom) contains two Cooperons, 
describing propagation in the intervals {t2, —ti) and {—ti,ti). 
Curved lines represent interaction propagators. 



where W^'i(ri4) = | /(dg)£|'e*«''". This is ultraviolet 
convergent only for d = 1, to which we henceforth re- 
strict our attention, but then W^^^ru) is infrared diver- 
gent. However, when Eq. (jBip is inserted into Eq. (|20|l . 
this divergence can be arranged to cancel between the 
two terms of Eq. I|B1(I : its second term, which stems 
from S*''''^, produces [using (7(7-14; ti, ti) — S{ri4)] a con- 
tribution — Wp'[(0) C(ri2; ti, i2), which can be rewritten 
as -W^'i(O) / driC {r 14; ti,-ti)C{r 42; -ti,t2), so that it 
takes a form similar to that resulting from the first (ver- 
tex) term of Eq. (jBl|l . Thus, the Bethe-Salpeter equation 
()2()|l can be written as 

{-DVl^ + dt, + 7ff) Ciri2;ti,t2) = 5(ri2) Situ) + 
dr4Wci{n4) C{ri4; ti,-h)C{r42; -h, t2) , (B2) 



where the kernel Wci(ri4) = W^'j(ri4) — Wci(*^) ^'^^ 
infrared problems: 



Wci(r) = 



dq Ae^T e^?'' - 1 
27r hui (p 



(B3) 



(Here ai = a^2ve?D is the inverse resistance per length 
of a quasi l-dimcnsional wire of cross section .) 
Eq. I|B2|) is solved by the following path integral: 



, AD 



Ciri2;h,t2) ^0{ti-t2) 



Vr exp 



-initi - 12) 



+ dt'W4r{t') ~ r{-t')) 
Jo 



(B4) 



where i* = ti for -t2 > ti > 0,tl = -t2 for ti > -t2 > 
and t\ = Q otherwise. Indeed, upon differentiating with 
respect to ti , we find that the contribution from the Wc\- 
terni in the path-integral is nonzero only for < ii < 
— as is the case for the IVci-term in Eq. (|B2|) . and in 



fact precisely equals the latter: 

VrW.Mti) - r{-h))eM- ■ ] 

dr4W,x{ri4)C{ri4]tu-ti)C{r42\ -tiM) -(65) 

The validity of this equality would be obvious if we were 
treating unmodified Cooperons. The fact that it remains 
true even in the presence of noise is due to the special 
nature of this noise: In principle, we need a "three-point 
Cooperon" involving times ti, —ti, and t2, and it factor- 
izes in the manner shown in Eq. (|B5|I only because the 
classical Nyquist noise correlator is a (5-function in time: 
indeed, as can be seen from Eq. ljB2p itself, the noise 
correlator always only connects time-points t' and — f', 
which means there is no correlator connecting two points 
in the disjoint intervals (<27~^i) and (— ii,ii) involved 
here (see Fig.O bottom). 

By rewriting the path-integral (|B4|I for the special case 
of equal times t2 = —rj, ti = 77, it can be shown to 
be identical to the exact Cooperon path-integral expres- 
sion considered by Altshuler, Aronov and Khmelnitskii 
in their seminal work (Eq. (20) in Ref. |3, correctedi^ for 
factors of 2), where (7(0; ?/, —77) — Tc\C^^{r, r), 



dt'W,Mt')-ri-t')) 



(B6) 



2e2T n r+°° dk 



{l-cos(fc[r(0-r(-0])} 



In comparing the expressions, note that our plays the 
role of their as a given extrinsic decoherence rate, 
and our the role of their r^r. 

Finally, we comment on the connection between our 
Bethe-Salpeter equation ljB2|l . which is quadratic in 
the Cooperon, and AAK's differential equation for the 
Cooperon (Eq. (23) in Ref. H correctedi^ for factors of 
2), which is linear: 



077 op^ 



fuTl 



S{V)S{P-P') 



V2- 



(B7) 



Here p{t) = [r{t) — r(— 1)]/\/2 is the difference coordi- 
nate introduced by AAK. Both equations yield the same 
result, since they are solved by the same exact path inte- 
gral. The origin of the difference between the two equa- 
tions is that Eqs. and (jBTj) is that AAK's Eq. ||B7|) 
describes the symmetric evolution of both the first and 
second time arguments of their Cooperon, namely 77 and 
—77 (as illustrated in Fig.|21 top), wheras in our Eq. (|B2|I . 
only the first time argument of the Cooperon, ti, is time- 
evolved. For 77 = ti = —t2, an integral such as the one 
in the second line of our equation (|B2|I becomes linear 
in (7, since the second Cooperon reduces to a S{r24 func- 
tion. However, such a simplification is possible only if 
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the interaction propagator is a i5-function in time, and 
therefore cannot be employed to simplify evaluation of 
the full Bethe-Salpeter equation in general, with interac- 
tion propagators more long-ranged in time. 



Appendix C: SOME SUBLEADING 
CORRECTIONS TO AND C'^'^ 



in Eq. HC3|) . which will be needed for d = 2, 3, and to can 
be taken as the elastic scattering time t^i. The fact that 
a cutoff is needed is of no great concern, since the lead- 
ing long-time behavior will turn out to be subdominant 
anyway. 

In the limit Tt ^ 1, we need only the asymptotic 
small-x behavior of Vd{x) in Eq. (|C3|) . which is given 
by 



In this appendix we calculate some subleading cor- 
rections to the decay function and the first order 
Cooperon C^^^'^ , which were mentioned but not discussed 
in detail in the main text. 



1. Calculation of Ff 



In Eq. 



in Section IIII Bl we expanded the prop- 



agator Cf^^q{u)) in powers of w, and subsequently eval- 
uated only the corresponding lowest order contribution 
Fq to the decay function, arguing that small lo dominate 
in the long-time limit so that the higher terms are neg- 
ligible. Let us now check this explicitly by calculating 
the first correction, F^^\t) = {F^^^(t))^, starting from 
Eqs. H32|) . and using methods and notations analagous 
to those of Paper I, Sections IIII PI and IVI Al It will 
be found to be subleading, so wc shall only calculate 
its order of magnitude, without caring about numeri- 
cal prefactors. We begin by noting that Eq. H29|l yields 

i^Tin),qi^))e = -W^\cu)/i^ , where 



pp 



W^p('^)=8T(«ech2[(s-'^)/27^])e 
Writing its Fourier transform as 



(CI) 



(C2) 



we note that the function w^{z) is dimensionless, 
peaked around zero, with heigth, width and weight all 
of order 1. Inserting this into Eq. H32|l for n = 1, writing 
the latter in a form similar to Eq. H33() . and rewriting 
the time integrals in terms of the dimensionless sum and 
difference variables x — t34/t and z = t^^T = xTt, with 
X = t^/t, we obtain 



^i'^W^^y^ dzw^^i{z)V^^\z/Tt) , (C3) 



where we defined the operator Dt = {—dt+d/{2t)+'yH), 
the function V^^\x) = J^^ ^dx 5p'^\t,t), with 



2t 



(dq) |i34|)-P('o7)(g- |i34|) 



T--d/2 



-d/2 



(C4) 



and the shorthand notations r = (1 — x)x, r = (1 — x)x. 
Moreover, we introduced an ultraviolet cutoff xq = to/t 



Vi'\x) 
V^'\x) 
H'\x) 



1 

1 



-1/2 



ln(a::a;o) 



-3/2 



2x, 



-1/2 



(C5) 
(C6) 
(C7) 



Inserting this into Eq. (|C3|) . we find that F'^\t) is 
smaller than the leading term F^^{t) =— (F^^j (t))^, 
given by Eqs. fII73() of paper I, by powers of the pa- 
rameters 7t = XjTt and 7// = which for present 
purposes are both 1: 



F^\r) 



F^'\t) = Fi^{t)0 ^H,lu{lHlult)/{Tto) 



FfP{t)0 

FFit) 
In(Tt) 



~ -1/2 -3/2 



1/2 



(C8) 

(C9) 
(CIO) 



In particular, recalling the leading time dependence of 
FfP(i) [namely t^/"^, tln{Tt) or t for d = 1,2,3], we see 
that the leading terms of Fj^''(t) all either vanish in 
the limit of no magnetic field (7// = 0), or are constant 
or decreasing functions of time. Hence we conclude 
that F^^\t) indeed can be neglected for the purposes of 
determining the decoherence time. 



2. Long-time behaviour of C^lifi 

Next, we consider in more detail the correction 
C^cifzit) = ^ici''/z(0; 5*7 to the Cooperon arising 

from 5]scif,z. j-j- jg giyen by an equation similar to the first 
term of l(?7jl . but using [Eq. ljA5d|l ] as self-energy. 

This contribution was purposefully omitted in our cal- 
culation of C^J^^'it) = C'£)^^(0;ii,-it) from (gH). To 
justify this omission, we shall now show that the long- 
time behaviour of C^l^fzit) is subdominant as compared 

to the leading behaviour from {cj^J^ (t)) ^. 

It turns out that for this calculation, we have to in- 
troduce both an IR cutoff in q, and an UV cutoff in 
uj [though the latter would not be necessary if, in con- 
trast to the calculation below, one would use the gen- 
eral expression (jA4p for the interaction propagator in- 
stead of the unitary limit l|A3b|l ]. The occurence of these 
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divergencies is not a surprise, since the self-energy di- 
agrams used in our Bethe-Salpeter equation constitute 
only a subset of the diagrams that make up the cross- 
terms of interaction- and weak-localization corrections to 
the conductivity, namely that subset of diagrams capa- 
ble of being iterated in a diagrammatic equation for the 
Cooperon. However, it has been shown in Ref. ^ that 
when all contributions to the conductivity to first order 
in the interaction (and second order in 1/g) are calcu- 
lated, numerous additional terms arise which turn out to 



cancel the abovementioned IR and UV divergencies, but 
which we have not considered here. 

For present purposes, it is sufficient to simply cut off 
these divergence by hand, since we shall find that the 
leading long-time behavior of this term is subdominant 
anyway. To identify this long-time behavior, we shall 
isolate the strongest singularity in the frequency domain 

of the Fourier transform C^^lllzi^) = / dte"^^C^^^Ilzit)' 
which has the following form: 



tanh(- 



2T 



^-tanh(£±|^) 



2T 



I[q,uj,uj] 



I[q,oj,u;] = / (dq)C°M2 [im^f + im^)]'') -C;%(c. + c^) 



(Clla) 
(Cllb) 

We shall analyse these expressions explicitly only for d — I, which is most prone to infrared divergencies for a; ^ 0, 
which would correspond to a long-time behavior that grows with time t; if C^l\{zi^) found to be subleading for 
d = 1, the same will be true for d = 2, 3. Using Eq. ljA3b|l for £.^{lj), the contour integral over q yields / = (/i — /2)/a^ 
[referring to the first and second terms in square brackets] , with 



1 



27 - i{2uj + u;) + 



i 



VD{2z2)[i^q + - 4? \/^(2zi)2[(zi - y^q)^ - z|] 



1 

Zl 



2(zi - ^/Dq) 



(zi - VDq)^ - zl 



(C12a) 
(C12b) 



where z\ — i^J^n — «w and Z2 — i\/^H — + uj). 

For 7/f — > 0^, the part of / which will yield the most 
singular contribution in the frequency domain (after in- 
tegration over uj and q) is a I/a/oJ singularity at w = 
(/i 2 both yield a a;"^/^ contribution, but those cancel in 
h'-h): 



1 



a'^{uj + iDq^f V0+ - iuo 



(C13) 



The occurence of a cj^^/^ singularity can be understood 
as follows: for uj — Q, the integrand in /[g, w, 0] diverges 
as l/{Dq^) (for d = 1), but this divergence is cut off by 
w 7^ 0, so that the integral goes as I[q,u}^uj] ^ u:'^!"^. 
(By a similar argument, it follows that for d = 2, 3, the 
leading w dependence will be less singular, namely \vl{lo) 
or a constant, respectively). 

The subsequent integrals of this term over q and u) 
need an IR- and UV-cutoff qm\^ and (iJmax, respectively. 



Keeping only the contribution that dominates for g,, 
and Wmax ^ 00, we find: 



(l)e 
selfZ 



In 



<e,T) 



ha'^i^Dy/Q+ - iuj VDqn,i 
The corresponding temporal behaviour is 



ln( 



max(e,T) 



Dq„ 



(C14) 



(C15) 



which is subleading at long times compared to cj^J^ (t) ~ 
[from Eqs. (gH) and HH)]. Even upon inserting 
gmin (Dt)-'^/^ for the IR cutoff, so that C'^J^lzi''^) ~ 
l/{ha^vD), we note that as long as gi(L^) ^ 1, this 
contribution is much smaller than the leading contribu- 
tion to (7(^^^'^(f) at scales t ~ r^, namely C'dec^('^v) — 
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